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Abstract 

We argue that the way to get the general solution of a Schrodinger equation in the presence of a time- 
dependent linear potential based on the Lewis-Riesenfeld framework is to use a Hermitian linear invariant 
operator. We demonstrate that the linear invariant proposed in p and q is an Hermitian operator which has 
the Gaussian wave packet as its eigenfunction. 

PACS: 03.65.Ge, 03.65.Fd 

The time evolution of a quantum system subject to a spatially uniform, time-dependent force has attracted 
considerable interest lately. The exact propagator for this system has long been known [1], as have a set of exact 
solutions (the Volkov solutions) [2-3]. Recent work [4-13] focuses further on exact solutions and their proper- 
ties. First, Guedes [5] obtained, by means of the invariant operator method introduced by Lewis and Riesenfeld 
(L-R) [14], a special solution for the time-dependent linear potential. The idea is that any operator satisfy- 
ing the quantum Liouvillc-von Neumann equation provides its eigcnstate as a solution of the time-dependent 
Schrodinger equation up to a time dependent phase factor. Later on, Feng [6] followed a method based on 
spatiotemporal transformations of the Schrodinger equation to get the plane-wave-type and the Airy-packet 
solutions where the one in Ref. [5] constitutes a particular case which corresponds to the so-called "standing" 
particle case in a linear potential [6] . 

However, Bekkar et al. [8] pointed out that the Airy-packet solution is in fact only a superposition of 
the plane- wave- type solution. In his Comment [7], Bauer explained that the solution found by Guedes [5] is 
simply a special case of the Volkov solution, with a zero wave vector fc, to the time dependent Schrodinger 
equation describing a nonrelativistic charged particle moving in an electromagnetic field. Dunkel and Trigger 
[10] considered the initial minimum-uncertainty Gaussian for a sinusoidally time-dependent linear potentiaL 
Bowman [11] investigated the time evolution of the general quantum state for the time-dependent linear system, 
which was shown to be that of a free-particle state, plus an overall motion arising from the classical force . Sang 
Pyo Kim [12] has shown that a charged free particle in a constant and/or oscillating electric field has a bounded 
Gaussian wave packet, which is a coherent state of a one-parameter and time-dependent ground state for 
the free-particle Hamiltonian. From the test-function method, Gengbiao Lu et al.[13] constructed an exact n 
wave-packet-train (n GWPT) solution , whose center moves acceleratively along the corresponding classical 
trajectory. 

Very recently, Luan et al. [9] have reexamined the linear invariant proposed by Guedes [5] and indicated 
that besides the solutions described in Refs. [5-8], a Gaussian wave packet eigenfunction solution can naturally 
be derived from the LR method if a non-Hermitian linear invariant is used, and they argue that this solution 
was ruled out before, because the authors in Refs. [ 5-8] assumed in advance the linear L-R invariant I{t) = 
A(t)p + B(t)x + C{t) be a Hermitian operator. 

The authors [9] stated that in Ref. [5, 8] , the Gaussian wave packet eigenfunction solution was ruled out , 
because the linear LR invariant is supposed to be Hermitian operator, and the B parameter was set to zero for 
hermiticity reason. In fact this is not the reason for setting B — in Ref. [8], it was meant to justify the 
B = imposed in Ref. [5]. We will show and we will correct in the present paper, that if B is different from 
zero {B ^ 0), we still obtain Gaussian wave packet (GWP) solutions. In fact, we could find all the results using 
the Hermitian invariant operator linear in p and q with any conditions imposed to get physically acceptable 
solutions as was done in Luan et al. paper [9]. In the latter approach, the authors obtained solutions labelled 
by complex parameters and complex eigenvalues, this is ambiguous. 



We recall that according to the theory of Lewis and Ricscnfcld [14], an invariant is an operator that must 
necessarily satisfy three requirements: 1. It is Hermitian. 2. It satisfies the Von Neumann equation. 3. 
Its eigenvalues are real and time-independent (we signal that the Hermiticity of the invariant is one of the 
essential conditions which makes the eigenvalues of the invariant time- independent). Furthermore, any invariant 
satisfying these three requirements leads to a complete set of solutions of the corresponding Schrodinger equation. 
So, a conventional solution is constructed as a linear combination of these solutions. 



The problem is to find the solutions of the Schrodinger equation, 



for the Hamiltonian 



d 

in—tp{x,t) = H ip{x,t) 



H{x,p,t) = —p^-F{t)x, 



(1) 



(2) 



where F{t) is a time-dependent function. 

According to the theory of Lewis and Riesenfeld [14], a solution of the Schro"dinger equation with a time- 
dependent Hamiltonian is easily found if a nontrivial Hermitian operator I{t) exists and satisfies the invariant 
equation 



dl dl 1 „ 



(3) 



Indeed, this equation is equivalent to saying that if 'f\{x, t) is an eigenfunction of I{t) with a time-independent 
eigenvalue A, we can find a solution of the Schrodinger equation in the form il>x{x, t) = exp [iax (t)] ^Pxix, t) where 
ax {t) satisfies the eigenvalue equation for the Schrodinger operator, 



d 



Tiax {t) ipx = i 

It turns out that the time-dependent invariant operator takes the linear form [9] 

I{t) = A{t)p + B{t)x + C{t). 



(4) 



(5) 



The invariant equation is satisfied if the time-dependent coefficients are such that 

B{t) = Bo, 



A{t)=Ao-^t, 
m 



C (t) =Co-Ao I F{T)dT + I F{T)TdT, 



(6) 
(7) 



lo / F{T)dT +^ [ 

Jo m Jq 

where Bq, Cq are arbitrary real constants. 

The eigenstates of I{t) corresponding to time-independent eigenvalues are the solutions of the equation 

Iit)^xi^,t) = >^ ^xix,t). (8) 
It is easy to see that the solutions of Eq. (8) are of the form 

'2{X-C{t))x- Bqx'^'^ 



^\{x,t) = exp 



2A {t) 



(9) 



Substituting Eq. (9) into Eq. (4) and accomplishing the integration, we obtain 
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«a(0)-/ {'-^^^^\dr-^Ln[J^], (10) 



[ 2mhA{Ty 2mA (r) 

\ 2mhA{Tf J \\l A{t)J' 

note that the logarithmic term goes downstairs as a time-dependent normalization factor in ip)^{x,t). 

Therefore the physical orthogonal wave functions ip^{x,t) solutions of the Schrodinger equation (1 ) are 
given by 



2iX-C(t))x- Bqx^ 
2 A (t) 



(11) 



A{t) ^ [ Jo 2hmA {t'f 
It is easy to veriiy that 

(^^aI^a') = ^(A-A'). (12) 
Furethermore, the evolution of the general Schrodinger state can be writen as 

/ + 00 
g{\)il^^{x,t)dX, (13) 
-oo 

where g (A) is a weight function which determines the state of the system such that {x, t) is square integrable 
i.e. J \ {x,t) \ dx is time-independent finite constant. 
In Ref. [8] , it has been shown that the choice 

5(A) = expf^') (14) 



leads to Airy function solutions. 

Here we shall show that one obtains a general Gaussian wave-packet solution by choosing weight function 

as a gaussian too, 

where a is a positive real constant. 

Substituting Eq. (11) and (15) into Eq. (13) and accomplishing the integration, we obtain the normalized 
Gaussian solution as 



*(x,i) 



\^%A{t)^(a + iSl^i-^dt'] 
-2C {t)x-BQx'^' 



■ exp < —I 



' Cjt'f 
2hmA{t') 



,dt' 



exp 



2A {t) 



exp ■ 



Jo RmA(t')^ ^ hAit) 



2 N 



(16) 



Let us now evaluate the mean value of x, p and the quantum coordinate and momentum fiuctuation in the 
state ^ {x, t) . After some algebra we find that 

{x) = {t) \x\ * {t)) = -A it) f ^^dt', (17) 

Jo mA{t') 

which is nothing but the classical position Xc{t), and 
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Ait) Jo mA{t'f 



is a classical momentum Pc{t). The position uncertainty 



Ax = ^(x2) - {xf = TiA it) \^ 
and the momentum uncertainty 



+ (/o IhmAit'f ^*') 



(18) 



(19) 



1 
Ax 



\ 



1 ^^,_n^BoAiti 



4a mA{t') 



(i) 2?ItoA (t')^ '^^ ) 



leads to the uncertainty relation 
ApAx 



\ 
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4a (t') 



2;imA(i')' 




> 



(20) 



(21) 



We can rewrite (16) as follows 



^{x,t) 



\ 



UA 



[t] {a 



■ * Jo 2hmA{t'y 



dt 



') 



exp<;- 



\/2^Aa;2 
-2C {t)x-BQx'^ 



exp < — i 



C(i 



7^2 



2;im^ (i') 



2A (i) 



exp 







rdt 



4aAa;2 



(22) 



Morever, the time-dependent probability density associated with this Gaussian wave packet is Gaussian for 
all times 



i*(^,i)r 



1 



27rAa; 



exp - 



{x-{x)f 
2Aa;2 



(23) 



we see that Aa; represents the width of the wave packet at time t. It is also readily verified tat the time-dependent 
probability density is conserved 



/oo 
{x,t)\^ dx = 1 
-OO 



(24) 



Equations (22) and (23) describe a Gaussian wave packet that is centered at a; = (x) whose width Ax{t) 

varies with time. So, during time t, the packet's center has moved from x = to x = —A (t) mMv)'^ ^^'^ 

1 

dt' ) . The wave packet 



its width has expanded from Axo = hAoy^ to Ax{t) = Axq^^^^J + ( /o 2nmaA(t')'^ ^ 
therefore undergoes a distortion; although it remains Gaussian, its width broadens with time whereas its height, 
, decreases with time. Further, it should be noted that the width of the Gaussian packet does not depend 



/27rAx 



on the external force F{t). Thus the shape of the wave packet is not changed by the external force. This means 
that the external force F{t) acts uniformly in the wave packet. 
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